The Propagation of Quantum Information Through a Spin System 
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It has been recently suggested that the dynamics ol a quantum spin system may provide a natural 
mechanism for transporting quantum information. We show that one dimensional rings of qubits 
with fixed (time-independent) interactions, constant around the ring, allow high fidelity communi- 
cation of quantum states. We show that the problem of maximising the fidelity of the quantum 
communication is related to a classical problem in fourier wave analysis. By making use of this 
observation we find that if both communicating parties have access to limited numbers of qubits in 
the ring (a fraction that vanishes in the limit of large rings) it is possible to make the communication 
arbitrarily good. 
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In quantum information science it is a crucial problem 
to develop techniques for communicating qubits. A key 
task is to develop protocols for taking a qubit's state at 
one location to another while minimising the degradation 
of the quantum coherence. 

In this paper we consider whether it is possible to use 
a one dimensional arrangement of qubits - a ring - cou- 
pled by nearest-neighbour interactions, to communicate 
a qubit from one part of the ring to another with high- 
fidelity. 

It is clear that this communication can be done, for 
example, if one can perform fast local unitary opera- 
tions on the individual qubits (by fast we mean that 
the local unitary operations are effectively instantaneous 
compared with the time it takes the interaction between 
qubits to change the state appreciably). This is because 
any nearest-neighbour interaction can be used to perform 
a swap operation on pairs of qubits J^, ^ JJ- Thus 
the communication can be performed by putting a given 
qubit into the state to be communicated and then mov- 
ing it along the chain of qubits by swapping the state 
along the line to the required position. 

In this paper, however, we have in mind the situation 
that we do not have access to the whole set of spins during 
the communication; rather the set of spins behaves more 
like a fibre into which one can input a qubit state, and 
hope that the state appears at the other end of the fibre. 
The question is whether it is possible for a set of spins 
with only its passive interaction to act as an effective 
communication channel. 

We would also like the qubit fibre to have the property 
that we could increase the length of the fibre without 
needing to completely re- fabricate it; thus if we wanted 
to communicate twice as far we would like to be able 
to simply add two identical fibres together and find that 
the fidelity and time of transmission depended in a simple 
way on the fidelity and time of the individual fibre. 

Some interesting work has already been done in this 
area concerning communication of qubits along a ID 



Heisenberg spin chain la (and, more recently, in the ID 
XY model, p). In [3 the aim was to put one qubit 
in a chain or ring into a given state and to transfer the 
state, as well as possible, along the chain (or ring) to 
the output qubit. The input and output qubits were 
placed at diametrically opposite places in the chain (or 
ring) and the interaction between neighbouring qubits 
was time-independent and constant around the length 
of the system. Interestingly if a ring has length four, 
there is a time after which the state of the output qubit 
is precisely the same as that of the input qubit. As the 
length of the system is increased, however, the maximum 
fidelity achievable was found to go down with increasing 
separation. For example, chains of length approximately 
80 have fidelity not much better that 2/3 (this is the fi- 
delity of the best classical transmission of an unknown 
qubit); thus the approach would not be useful for trans- 
mission over long distances. A further notable feature 
of the protocol is that the time at which the maximum 
fidelity is reached does not depend in any simple way on 
the length of the ring. This is because the maximum 
fidelity arises when an expression of the form e""^* 
reaches its maximum (where t is the time, and the aj 
are rational constants which depend on the length of the 
ring). Thus in this framework there is no natural notion 
of the speed of transmission along the fibre. 

Further progress was made in where the authors 
show that a single qubit can be communicated perfectly 
in a hypercube geometry. The authors also show how to 
communicate a qubit perfectly in a linear chain with vari- 
able couplings given by a specific arithmetic sequence. 
Using this idea for communicating arbitrary distances, 
would however require re-fabricating the channel for each 
distance. 

A further interesting proposal for communicating 
quantum information using the ground state of a certain 
spin system has been recently developed in 0|. There it 
was found that perfect transmission is possible when the 
freedom to apply local measurements on all the spins is 
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allowed. 

We can summarise our desiderata in this paper by 
saying that we wish to develop schemes for high fidelity 
quantum communication with three essential criteria; (i) 
Minimum control requirements. We require that it is un- 
necessary to apply many fast control operations through- 
out the transmission; (ii) Robustness. We want the de- 
vice to be capable of tolerating small errors without di- 
minishing the quality of the communication too much; 
and (iii) Flexibility. We want to be able to change param- 
eters (such as the locations of the communicating parties) 
simply without requiring a new device be fabricated. 

Our proposal makes quantum spin chains behave as a 
medium through which quantum information propagates. 
Using the fact that dynamics in restricted subspaces for 
ID quantum spin chains can be described by classical 
fourier wave analysis we show that there are natural no- 
tions of group velocity and dispersion for quantum infor- 
mation. We propose a figure of merit for the information 
propagation, the qubit rate, and analyse this quantity for 
a family of ideal and realistic physical models: nearest- 
neighbour ID spin chains. In addition, we argue that our 
protocols satisfy the requirements (i)-(iii) we outlined in 
the previous paragraph. 

We show how to increase the qiiality of quantum com- 
munication in a variety of realistic spin models by relax- 
ing the restriction that the two parties can only access a 
single qubit. We allow both parties the ability to access 
a limited number of extra sites, a fraction which tends 
to zero as the length of the chain increases. The con- 
nection between propagation of quantum information in 
spin rings and classical fourier wave analysis provides a 
way to visualise the propagation of quantum information 
pulse through a medium. 

We now turn to a precise description of the protocol 
we will study. We consider a ring of N spin-1/2 systems 
evolving according to a nearest-neighbour hamiltonian 
H. We identify the {N + l)th site with the 1st site, i.e. 
An+i = Ai. We imagine that two parties, Alice A and 
Bob B are located at sites 1 and N/2, respectively. (We 
assume, for simplicity, that the ring is composed of an 
even number of spins. Note, however, that our subse- 
quent results do not depend on this fact.) We suppose 
that Alice and Bob are able to access a total of A sites 
each, centred on the 1st and N/2th sites, respectively. 
Alice and Bob are allowed to perform any operation al- 
lowed by the rules of quantum mechanics on their sites. 
Finally, we assume that the completely polarised state 
|00 • • • 0) is an eigenstate of the hamiltonian H, and that 
before the protocol begins the system is prepared in this 
state. 

Alice wants to communicate a (possibly unknown) 
qubit state \tp) to Bob. Given a specification of H (as- 
sumed to be known to both Alice and Bob) Alice per- 
forms some encoding operation Ua on her A qubits. The 
system is now in the state |^'(0)) = Ua\00- ■ -0). As the 



ring is always evolving according to H this state imme- 
diately begins to evolve: \^{t)) = e"'-^*|*(0)). Bob now 
waits a certain period of time T which depends only on 
H, A and N, and when this duration elapses he performs 
a decoding operation Ub on his addressable qubits in or- 
der to decode or refocus the communicated state into one 
of the qubits in the ring. (Additionally, Bob could ap- 
ply a swap operation to move the decoded state into a 
static register qubit. We prefer to ignore such register 
qubits because we imagine that decoding unitary oper- 
ations could also be performed as intermediate steps in 
a kind of quantum repeater.) The protocol is deemed to 
succeed when the average fidelity 

T{H,N,K,T,Ua,Ub) = ^J dn{ij\p\ij), (1) 

where \ip) is the input state, which we average over the 
Bloch sphere, and p is the decoded output state, is above 
some prcspccified threshold value t. 

Ultimately we are interested in how well a spin ring 
performs as a conduit for quantum information. In this 
more general scenario we allow Alice and Bob the free- 
dom to apply a number of encoding/decoding operations 
in succession. The objective is to maximise the number 
of qubits successfully communicated (i.e. when > t) 
per unit time. Write the maximum number of qubits that 
can successfully communicated (where the maximum is 
taken over the encoding/decoding operations) in a du- 
ration T as Mt{T). We define the long-time average 
Qr{H) = limT^oo t^Mt{T) to be the qubit rate for H. 

Obviously the evaluation of Qt-{H) is extremely diffi- 
cult, even for the simplest systems. We only obtain lower 
bounds for this quantity for a class of rotationally invari- 
ant nearest-neighbour hamiltonians. 

There is a simple case where we can evaluate Qt{H) 
exactly, namely when if is a hamiltonian corresponding 
to the translation operator H = — ilog(T), where T is 
defined by the following action on computational basis 
vectors, T(|ai, 02, . . . , un)) = \a2, as, ... , ajsr, oi). When 
H is of this form the qubit rate takes the maximum value 
Qr{H) = 1 for all r. (We have adopted rescaled time 
units; we will discuss how to calculate the constant of 
rescaling later.) In order to see this first note that when 
t = n is an integer the propagator U{n) = e"*-^" = T" is 
a power of the translation operator. To achieve the qubit 
rate, Alice needs to encode the state of a qubit into one 
site of the ring at every integral t. 

As we'll argue presently, every hamiltonian H corre- 
sponding to the translation operation is massively non- 
local and contains interaction terms between many sep- 
arated subsystems. For this reason it is unlikely that a 
system will be fabricated which naturally evolves accord- 
ing to H. 

Consider now the translation operator T. It is easy to 

see that the action of T on computational basis states 
breaks up into blocks, for example, the states |00---0) 



3 



and form blocks all by themselves, and the 

states j ^ 1, . . . ,N, where |j) denotes the state of all 
|0)'s with a |1) in the jth site, form a block (we'll often 
refer to these states as one-particle states). The remain- 
ing basis states each segregate into blocks in a similar 
fashion. Consider a block of states of size M. Choose a 
state from this block and call it \ai). Every state \ak) 
within the block can be written as \ak) — T'^~^|q!o). Us- 
ing these M states one can construct M eigenstates of 
T: |/3fe) = l/VMY^f^^v^^-'^^^\aj), where v = e^' is 
the Mth root of unity. We note that each set of eigen- 
states so constructed is a quantum fourier transform of 
the set {|aj) I j = 1, . . . , M}. Performing this procedure 
for each block gives rise to the complete set of eigenstates 

for r. 

We study the qubit rate for the class of rotationally 
invariant nearest-neighbour spin-1/2 models. Because of 
the rotational invariance [H, T] = and we may simul- 
taneously diagonalise both H and T. We restrict our 
attention further to the class of ID spin rings which fix 
the following special eigenstates of T, which we call the 
twisted W-states, \W{k)) = ^ Ej^i M^^^^^'^b')' where 
/i is the A^th root of unity /i — e^'\ All the protocols we 
describe in this paper take place in the subspace formed 
by the twisted W^-states. 

The spin hamiltonian H fixes the states with 
eigenvalues uj{k) (and also, via a rescaling of energy, 
H fixes |00- • -0) with eigenvalue 0). Suppose we start 
the system in the arbitrary state |^(0)) — ajOO- • - 0) + 



/3E 



N 



N 



__,cj\j), where \a\' + = 1 and ^^-L^ |c,f = 1. 
The subsequent time evolution of this state can be writ- 
ten 



N 

|*(t)) =a|00•••0)-H/3^Cfce- 
fc=l 



I l¥(fc)), (2) 



where 



1 

the |j) term can be found as 



(j-i)'=c^. The coefficient Cj{t) of 



1 



N 
fc=l 

Because the dynamics only take place in the zero- and 
one-particle subspace we can establish the following use- 
ful result. Suppose a = 0. If we bipartition the spin sys- 
tem into two subsystems A and A we can always write 
the system's state \^{t)) 
way. 



S,=i 6^(0 li) in the following 



\^{t)) = ^ 1 - C-^{t) \<ty)A\Q)A + ^C^)\ 0) A 1 0') A> (3) 
where 



and 



10') 



Note that Eq. Q is a two-term Schmidt decomposition 
of \^{t)) because a{(I)\0)a = a{<P'\^)a = 0- 

When Alice wants to send the state \tp) ~ a\Q) + /3|1) 
she will encode this state as a|00 • • • + PY^^=i Cj(0)|j), 
where she is free to choose nonzero Cj(0) as long as the 
index j lies within her subset of addressable spins. Af- 
ter time T has elapsed then this state can be written, 
following the discussion in the previous paragraph, as 



\^{T)) ^ P^l ^Cb{TMb\Q)b+ 



|0)^(a|0) 



/5v/CMT)k/)B), (4) 



where now the bipartition is given by BB^ where B is 
Bob's addressable spins and B now refers to all the other 
spins and 



Bob now applies a decoding unitary Ub to his part B of 
the ring. The unitary is (partially) defined by Ub\0)b = 
I Ob) and UB\r]')B = \N/2). (There is a great deal of 
arbitrariness in how Bob decodes his state. We choose to 
concentrate the state into one qubit in order to facilitate 
the evaluation of the average fidelity of the channel.) 

After decoding, the state p of the qubit N/2 has density 
operator 



'12(1 -Cs(T)) y/CB(T)a(3' 
Wl'CBiT) 



(5) 



with respect to the basis {\Q)n/2j 11)^/2} of qubit 
This state may be written p — £t{\'4') n /2{''P\) = 
Mo|V')iv/2(V^|M(I + Mi|7/;)jv/2(V'|M/, where |^A)a,/2 = 
o-\^)n/2 + /3|l)Ar/2 and Mq and Mi are the Kraus op- 
erators of an amplitude damping channel £, 



Mo 



1 



and 



Ml 







^1-Cb{T) 



with respect to the basis {|0)jv/2, |l)jv/2} (cf- p)- 

The fidelity {ip\UBpUg\tp) of Bob's state with the input 
state is given by 

\a\^ + {1 + 2^C^)-CB{T))\a\^m' + CBiTM^ (6) 
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The average fidehty Eq. Q can be evaluated as | + 
i y/CB{T) + \Cb (T) . This function depends monoton- 
ically on the quantity Cb{T). 

It is clear from the discussion in the preceding para- 
graphs that the quantity Cb{T) plays a central role in 
determining the effectiveness of the communication pro- 
tocol. We now study this quantity further, and show 
how to reduce the problem of its maximisation to a well- 
known type of problem in Fourier wave analysis. 

Motivated by the importance of C_b(T) we introduce a 
method for visualising quantum states which lie within 
the one-particle subspace. Given such a state = 
a|00 • • • 0) -f /3 J2f=i Cjb), we visualise it by plotting the 
quantities i^j == |cjp against site number. Clearly the 
area under the resulting (bar) graph equals 1. Obviously 
this pictorial representation does not represent any of 
the phase information contained in the system's state. 
However, this phase information plays no part in the 
quantity we are interested in, Cb (T) , which is simply the 
area under the graph of the state \i^{T)) between sites 
N/2 - A/2, . . . , N/2 + A/2. The graph of the system's 
state |i^(r)) is given by 

1 ^ 

fc=i 

It is worth emphasising that the method we have 
just introduced for visualising quantum states of the 
ring depends crucially on the properties of our spin sys- 
tems. The dynamics for the ring occurs solely within 
the zero- and one-particle subspaces spanned by |0) and 
{|i) \j — 1, • ■ • ,N}. This subspace scales linearly with 
the ring. It is precisely this feature which facilitates our 
construction for the graph of the state, which is essen- 
tially a method for representing, visually, a linear number 
(in N) of degrees of freedom. 

The problem of maximising the area under the graph 
yj{T) between certain sites is a well-known a classical 
problem in fourier wave analysis; our problem has re- 
duced to studying the linear dynamics of superpositions 
of complex scalar waves. The quantity u){k) is the dis- 
persion relation and the position variable is given by 
Xj = {j — l)/N. We wish to work out an assignment 
Cj(0) of initial amplitudes in Alice's subsystem, or a 
wavepacket of quantum information, so that subsequent 
evolution preserves as much of the width and integrity of 
the wavepacket as possible. 

Given a specification of the dispersion relation w(fc) 
we have all the information we require to embark on the 
design of wavepackets which preserve their shape. These 
numbers are, of course, the eigenvalues of H for the eigen- 
states \W{k)). 

To illustrate our ideas in the following we use the 



Heisenberg model on N sites 



where we have introduced an arbitrary constant x- We 
will explain which value of x we choose in the follow- 
ing. The dispersion relation for the Heisenberg model is 
u;(fc) = 2x (l — cos(2^)). (We have rescaled the zero of 
energy for the Heisenberg model so that |0) has eigen- 
value zero.) 

Currently our problem is not identical to the well- 
studied problems in wave motion (see, for example, 9]). 
The difference is that we are dealing with wave mo- 
tion in a discrete system whereas most of the theory of 
wave motion (at least in regards to the design of optimal 
wavepackets) is for continuous systems. In the subse- 
quent discussion we will make use of the results from 
the continuous theory, but it must be understood that 
these results only hold asymptotically in the limit where 
N oo and the intersite spacing S —^ 0. In particular, 
we emphasise that the appearance of derivatives must 
be understood as a formal device and really only make 
sense in this limit. We'll justify the applicability of these 
asymptotic results presently by comparing them with nu- 
merical results. 

In classical fourier wave theory the concept of group 
velocity plays a crucial role in the design of wavepack- 
ets which preserve their shape as time evolves 0- The 
group velocity v{kQ) for a wavenumber feg is defined to be 

^(^o) = ^ ^'^dk^ I k=ko • ^^^^ models we consider we 
are given a specification of uj{k) which makes sense for 
fractional k, so that it is possible to compute this deriva- 
tive.) The significance of this quantity is that an ini- 
tial wavepacket J2k=i WkC^'^''^^^ consisting of fundamen- 
tal waves with large amplitude focussed on wavenumber 
/co will translate with velocity w(/co) as time evolves. 

Unless the system is dispersionless (w(fc) = ck) the 
initial wavepacket will spread as time evolves Q. The 
rate-of-spread of the initial wavepacket is proportional to 
the second order term 77- ^S^^ I , , . If this derivative is 

271- dk-' \k=ka 

zero (as will be the case in one of the models we consider) 
there is no appreciable spread of the wavepacket to sec- 
ond order in uj{k). There may, however, be contributions 
from third and higher order terms, which cause an initial 
wavepacket to spread. 

The generic strategy 10] for designing wavepackets 
which change shape as little as possible is to use a 
gaussian-modulated wave of variance A centred at po- 
sition Xfe and wavenumber fcg corresponding to the max- 
imum available group velocity i'(fco). In our case the 
"wave" is a twisted VF-state, so this "wavepacket" is writ- 
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ten |*(0)) = a|0) + [3\G{xk, ko, A)), where 

N 



i=i 



and J\f is chosen to normahse the state. The graph for 
the initial state is given by 



N 

E 



(x ■ -Xj)-^ 

g h2irifcoa:3 



This state has exactly the form of a gaussian-modulated 
complex scalar wave centred around wavenumber fco and 
with (spatial) width L. (We define width L to be when 
more than some prespecified area, say 95%, lies within 
an interval L centred on the maximum of the gaussian. 
We note that L is some multiple of A. Throughout the 
remainder of this paper we choose L = 4A.) 

It is known that an initial gaussian wavepacket, un- 
der time-evolution, remains, to goo d approximation, a 
gaussian of width L{t) satisfying 



goo d 

m 



m 
m 



Lo"{ko)t 
L2(0) 



(7) 



Because Alice is assumed able to address only NL{Q) 
sites (recall we have rescaled the length of the ring to 
be 1 by defining the position variable Xj = {j — l)/N) 
we propose using the gaussian initial pulse introduced 
in the previous paragraph truncated (and renormalised) 
after 7VL(0) sites. 

We want to design a wavepacket that spreads by an 
amount that scales favourably with N. We demand that 




2000 3000 
Number of sites 



5000 



FIG. 3: Dots: size of the final wavepacket in the Heisen- 
berg spin ring for N = 50 to N — 5000. (Initial wavepacket 
is a gaussian modulated wave of variance Ns truncated af- 
ter 2 A*' 3 sites.) Final width F is defined to be when more 
than 95% of the wavepacket lies between F sites. Solid line: 
graph of 2.8N3 for comparison. The apparent jumps of inte- 
ger amounts greater than 1 in the numerical results is an arti- 
fact of the way the width of the final wavepacket is calculated. 
The final wavepacket has a sequence trailing oscillations. In 
order to get an area of more than some arbitrary amount (like 
0.95) it is sometimes necessary to include several more sites 
in a go. 
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a wavepacket of width L spreads by only a constant fac- 
tor to a final width of L' « kL by the time it reaches 
Bob. To simplify the discussion we note that for typical 
wavenumbers close to we have u;"(fco) = A/iV^ (this is 
true of all the models we consider) . We analyse the scal- 
ing, with t Ki N (the time it takes for a wavepacket to 

traverse the ring), of the spread S{N) of the wavepacket: 

1 

2 

1 ■ ' ■ 



S{N) = L{N)/L{0) 



X 



If we want the 



final width of the wavepacket to be independent of the 
length of the ring we need to choose i(0) = 1/^/N. This 
means that the initial packet needs to be supported on 
a subsystem containing on the order of ^/N sites. In the 
"mesoscopic" limit of large but finite N the fraction of 
the ring that Alice and Bob need to access can be made 
arbitrarily small (it scales as A^^a). 

We now analyse these results for the Heisenberg model. 
We want a signal pulse to travel as fast as possible. Look- 
ing at the group velocity uj'{k) — ^sin(^fc) we find 
that this occurs for k = 3A^/4, iy{k) = At 
this point we choose x — 1/4. We do this so that the 
final wavepacket will arrive at Bob's portion of the ring 
at time t = N (recall that Alice and Bob are separated 
by a distance N/2). So we choose the initial wavepacket 
to be a truncated gaussian-modulated superposition of 
width \/iV of twisted M^-states centred around k — N/i. 

The case of wavepackets centred on fco — N /A is rather 



Fig. (0) we see that the shaped pulse retains its shape for 
much longer durations. This is because fewer waves are 
involved in the superposition. We note that, generically, 
for wavepackets centred around wavenumbers not equal 
to A^/4 a square-root scaling for the final wavepacket is 
observed. 

Before we conclude our discussion of the numerical re- 
sults we refer to Fig. ©. In this graph we have plotted, 
for the Heisenberg ring, the width of the final wavepacket 
(we define the width F to be when more than 95% of 
the area lies within F/2 sites either side of site N/2) for 
rings of sizes 50 through 5000 sites. We have also plotted 
2.8A^3 for comparison. Because the Heisenberg model 
is diffusionless to second order in the dispersion relation 
Lu"{ko) — 0, we expect that the numerically recorded 
spread ought to be smaller than the spread predicted for 
models with second-order dispersion. We can see, from 
Fig. lO that this is indeed the case. For N = 5000 the 
proportion of the ring that Alice and Bob must be able 
to access is less than 1% each. 

To conclude we construct all the rotationally invari- 
ant nearest-neighbour spin hamiltonians which fix the 
|0) and twisted I^-states. An easy way to construct all 
such hamiltonians is to first assume that the hamiltonian 
preserves = X^jLi '^j ■ Any such hamiltonian may be 
written as 



iiitj case oi wavtjpcicKt'LS ceiibicu uii kq — iv/'i is launei jY 

special ~ note that the dispersion w"(fco) for this model ff = co/+^ \ci(7+ crj +C2CtJ (7+ +di{a+ aj^^+crj ct+^^) 



is identically 0. Naively applying formula Eq. (0 for 
the spread of the wavepacket suggests that an initial 
wavepacket of any width will not spread at all. This 
is an artifact of our approximation; we assumed that the 
dispersion relation could be expanded and truncated at 
second order with little error. In the special situation 
where oj'^ko) vanishes we need to go to third order. For 
this case a gaussian pulse does not remain a gaussian. 
(Instead, it becomes an Airy function The formula 

Eq. (|7|) is invalid in this case and we need to use the 
third- order broadening factor 



m 



i + - 



1 f iu'"{ko)t 

2 \V2L^iO) 



where = ^(cj ± icrj). It may be verified, with a 
little algebra, that this is the most general rotationally 
invariant nearest-neighbour hamiltonian which nontriv- 
ially fixes the twisted W^-states. (There is one additional 
term which fixes the twisted VF-states, "'^''^j+i' however, 
it trivially annihilates these states and thus contributes 
nothing to the dispersion relation. Note also that the /i 
term annihilates the twisted M^-states.) 

The hamiltonian Eq. ||SJ) gives rise to the following dis- 
persion relation 



We solve this equation for t = to find the width L(0) 
of initial wavepacket which has constant spread with N. 
We find that L(0) need only scale as N~3. Such a 
wavepacket consists of NL{Q) = N3 sites. 

We illustrate these results for the Heisenberg ring of 
A^ — 100 sites in Fig. (P) and Fig. ||2Jl. By way of con- 
trast, in Fig. (0 we first show the dynamics of the state 
with a single |1) at site N/2 (this type of localised initial 
state is what is used in the protocols of 0, and Q). 
Because the state \N/2) is an equal superposition of all 
the twisted VF-states all the waves are involved in the dy- 
namics and hence the initial packet disperses rapidly. In 



u}{k) = A + Bcoa 



27r. 
— / 

N 



S'sin ( ^fc 



(9) 



where A = cq + ci[N - I) + 02 + 2^2, B = 2di and 
B' = —2e\. We note that both the Heisenberg model in 
a magnetic field and the XY model in a magnetic field 
can both be expressed as in Eq. (jSJ for specific choices 
of Cj, dj, ei and /i. The dispersion relation Eq. (jSJ is no 
more general than that for the Heisenberg model so all of 
our discussion concerning the design of optimal signals for 
the Heisenberg model carries through straightforwardly 
for the general class of models Eq. © . 
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We see from the dispersion relation Eq. Q that 
the most general rotationally invariant nearest-neighbour 
hamiltonian fixing the twisted VF-states is dispersive. 
This shows us that, in order for a hamiltonian H to ex- 
ponentiate to the translation operator T = e'^, it is nec- 
essary for H to contain interaction terms between sepa- 
rated parties. This verifies our earlier claim that such a 
hamiltonian must be nonlocal. 

Finally, we note that our results imply that the qubit 
rate for a rotationally invariant nearest-neighbour Hamil- 
tonian on N sites is bounded below: Qo.9b{H) > 

We have introduced a method for improving the quan- 
tum communication characteristics of ID quantum spin 
rings. According to the connection between the dynam- 
ics of quantum information in these rings with fourier 
wave analysis we have been able to import many of the 
results concerning the design of signals which disperse 
minimally. Clearly our results illustrate that our commu- 
nication protocol has minimal control requirements and 
is fiexible. We have not tackled the problem of determin- 
ing the resistance of our protocol to error. 

Many future problems suggest themselves at this stage. 
Perhaps the most interesting is the extension of the re- 
sults to linear chains. In this case our results cannot be 
applied because the twisted VF-states are not eigenstates 
for nearest-neighbour hamiltonians on a chain. However, 
a generalisation of our method of pulse shaping ought to 
be possible. Another important future problem concerns 
the resistance of the protocol to errors. We believe that 
the protocol is robust, but a full analysis needs to be 
done. Finally, we note that all our protocols take place 
in the single-particle subspace. The dimension of this 
subspace increases linearly with increasing numbers of 
qubits; however the hilbert space of the system increases 



exponentially with N . Perhaps it is possible to increase 
the qubit rate substantially by taking advantage of larger 
subspaces which include two and higher particles. Fur- 
ther investigations along these lines are being conducted. 
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